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^ ■ Abstract 



We analyze the double radiative I?-meson decays B s — > 77 and i?^ — > 77 in 
QCD factorization based on the heavy-quark limit m& ^> Aqcd- We systemati- 
cally discuss the various contributions to these exclusive processes. The dominant 
effect arises from the magnetic-moment type transition b —>■ s(d) , ~y where an ad- 
ditional photon is emitted from the light quark (one-particle reducible diagram). 
The contributions from one-particle irreducible diagrams are power suppressed. 
We argue that they are still calculable within QCD factorization. They are used 
to compute the CP-asymmetry in B — > 77 and to estimate so-called long-distance 
contributions in B and D — > 77. Numerical results are presented for branching 
ratios and CP asymmetries. 



1 Introduction 



The rare decays of B mesons provide an excellent means to further explore the 
Standard Model (SM) and to detect New Physics. In particular the Cabibbo- 
favoured radiative b — > S7 modes belong to the small number of rare decays that 
are experimentally accessible already at present |]-{|] and give severe bounds on 
the parameter space of New Physics scenarios. On the theoretical side both the 
inclusive |5|-|7| and exclusive [jSHTOH b — > 57 decays are now known at next-to- 
leading-logarithmic (NLL) accuracy. 

For the double radiative decays B — > 77 experimentally so far only upper 
limits on the branching fractions exist: 

B(B S ^jj) < 1.48 ■ KT 4 at 90% C.L. [11] m 
B(B d -> 77) < 1.7 ■ 1(T 6 at 90% C.L. [pj ^ 

The Standard Model expectations are roughly two orders below the current upper 
limits. These decays have a rather clean experimental signature. Furthermore 
they are of interest because they could provide useful tests of QCD dynamics 
in B decays. The B — > 77 modes realize the exceptional situation of nontriv- 
ial QCD dynamics related to the decaying B, in conjunction with a completely 
nonhadronic final state and simple two-body kinematics. In principle they also 
probe the CKM parameters V ts and V t d and could allow us to study CP violating 
effects as the two-photon system can be in a CP-even or CP-odd state. While 
contributions of New Physics to B s — > 77 are strongly constrained by the mea- 
surements of b — > S7 transitions, a sizable enhancement of B d — > 77 from effects 
beyond the Standard Model is still conceivable. 

The theoretical treatments of B s — ► 77 performed so far (TT3|-[22j all had to 
employ hadronic models to describe the B s meson bound state dynamics. A 
clear separation of short- and long-distance dynamics and a distinction of model- 
dependent and model-independent features were therefore not possible. This 
concerns especially the dynamics of the light-quark constituent inside the B, but 
also contributions from intermediate J/ip, r) c , or meson states, which have 
been discussed as sources of potentially important long-distance effects. Another 
process related to B — > 77 is the decay B — ► Z + Z~7. This decay has recently been 



studied in |23| 



We present in this paper a systematic analysis of the exclusive double radiative 
decays B s d — ► 77 in QCD, based on the heavy quark limit m& 3> Aq C d- This 
limit allows us to give a factorization formula for the hadronic matrix elements 
of local operators in the weak Hamiltonian f8|,|9|, 10|,p4 1. In this manner we can 



systematically separate perturbatively calculable hard scattering kernels from 
the nonperturbative 5-meson light-cone distribution amplitude (LCDA). Power 
counting in Aq C D/ m b allows us to identify leading and subleading contributions to 
B — > 77. Only one diagram contributes at leading power, but the most important 
subleading contributions can also be calculated. The inclusion of these corrections 
is used to estimate CP asymmetries in B — » 77 and one-particle irreducible two- 
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photon emission from light-quark loops in D — > 77. These corrections represent 
the quark-level analogue of so-called long-distance contributions in these decays. 

The remainder of this paper is organized as follows. In section we give the 
basic formulas necessary for the calculation of the observables of interest. Section 
|3| contains our discussion of leading-power and section |] the one of subleading- 
power contributions. In section [5] we comment on long-distance contributions in 
B — ► 77 and D — > 77 decays. Numerical results are presented and discussed in 
section Section [7| contains our conclusions. 



2 Basic formulas 

The effective Hamiltonian for b — > 377 is identical to the one for b — > 57 transi- 
tions. This is because the equations of motion (EOM) can be used to reduce the 
basis with operators containing two photon fields ^5f . Up to corrections of order 
1/Myy the effective Hamiltonian thus reads 



H 



eff 



— V AW 
V2 ^ p 



p=u,c 



where 



The operators are given by 



C1Q1 + C2Q2 + 2! CiQi 



i=3,. 



+ h.c. 



X (s) 



V* s V p b 



(2) 



(3) 



Qi = 
QI = 
Qs = 


(sp) V - A (pb) V -A 

(SiPj) V -A(Pjbi)v-A 

(s6) y ^5I(gg)y-A 

q 


(4) 
(5) 
(6) 


Qa = 


q 


(7) 


Q 5 = 


(s6) y _A^(gg)y+A 

q 


(8) 


Q 6 = 


(Sibj)v-A^2(QjQi)v+A 

q 


(9) 


Qr = 




(10) 


Qs = 




(11) 



The most important operators are the magnetic penguin operator Q 7 and the four- 
quark operators Q\ 2 - The sign conventions for the electromagnetic and strong 
couplings correspond to the covariant derivative = <9 M + ieQfA^ + igT a A^. 
With these definitions the coefficients C 7j8 are negative in the Standard Model, 
which is the choice generally adopted in the literature. The effective Hamiltonian 
for b — ► g?77 is obtained from (PHI]]) by the replacement s —>■ d. The Wilson 
coefficients Cj in (|2|) are known at next-to- leading order 0. 
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The amplitude for the B — > 77 decay has the general structure 



= ~7=^-fB \A+ [2ki -e 2 k 2 -ei - m\t x • e 2 ) - 2z k 2 , e x , e 2 )l (12) 

Here and are the photon field strength tensor and its dual, where 

F>" = l -e^ Xp F Xp (13) 



with e 0123 = — 1. The decay rate is then given by 

r(5-77) = ^P^(K| 2 + |A_| 2 ) (14) 



In the heavy-quark limit we propose a factorization formula for the hadronic 
matrix elements of the operators in 

( 7 (e 1 )7(e 2 )|g^) = fs ['dCTnO^BiO^, (15) 

where the £j are the polarization 4- vectors of the photons and = is 
the leading twist light-cone distribution amplitude of the B meson. The latter 
quantity is a universal, nonperturbative object. It is defined via the light-cone 
projector for the B meson at leading power p4j 



(0\b(0)q(z)\B(p)) = ^(^+ m6 ) T5 [ 1 dte<&+*-[$ B1 (t)+ (16) 

4 jo 

The functions $bi,B2(0 describe the distribution of light-cone momentum frac- 
tion £ = l+/p+ of the spectator quark g with momentum Z inside the 5 meson 
(only $£i contributes in 1^). Here light-cone components of four- vectors t> are 
defined by 

V± = ^T (17) 

The wave functions are highly asymmetric with £ = 0(AQCD/ m b)- They are 
normalized as 

/V$m(0 = l /V<WO = (18) 
jo jo 

If the light-like vector n in fll6|) is chosen appropriately parallel to one of the 
4-momenta of the photons, only the first negative moment of $bi(£) appears, 
which we parametrize by a quantity Xb = 0(Aqcd), i-e. 

^•sM™ (19) 
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Because there are no hadrons in the final state, only one type of hard-scattering 
kernel T (type II or hard-spectator contribution 0J24J]) enters the factorization 
formula. The QCD factorization formula fll5|) holds up to corrections of relative 
order KqcD/ m b- The form of (|T5| ) with a simple convolution over the light-cone 
variable £ is appropriate for the lowest (leading logarithmic) order in a s , which we 
will use in the present analysis. A generalization to include transverse-momentum 
variables is likely to be necessary at higher orders in QCD p4] , p6| , ^7| . p8 | . 



We conclude this section with a brief discussion of CP violation in B — > 77. 
The subscripts ± on A±, defined in (|l"2"D for B — > 77, denote the CP properties of 
the corresponding two-photon final states, which are eigenstates of CP: Acp=+i 
is proportional to the parallel spin polarization e± ■ €2 and Acp=-i is proportional 
to the perpendicular spin polarization e\ x e* 2 of the photons. In addition to A± 
we introduce the CP conjugated amplitudes A± for the decay B — > 77 (decaying 
b anti-quark). Then the deviation of the ratios 

- H4+!4 (20) 

from zero is a measure of direct CP violation. We shall focus here on this direct 
effect, which is specific for the B — > 77 decay, and will not consider CP violation 
originating in B B mixing. 

Due to the unitarity of the CKM matrix we can parametrize 



A± = \ u a±e ia " + \ c afe ia < (21) 
A ± = A>±e^ + \* c afe ia < (22) 



where a^ c are real hadronic matrix elements of weak transition operators and the 
a^ c are their CP-conserving phases. We then have 

± = 2 Im(\ u \* c )a±at sm(a± - aj) , . 

CP |A u | 2 a± 2 + |A c | 2 a c ±2 + 2Re(A u A*)a^cos(«±-«J) 1 ) 

The weak phase differences read 

Im(A^A^) = T |A^||A^|sin7 (24) 

with the negative sign for q = s and the positive sign for q = d. The strong phase 
differences sin(aj — aj) arise from final state interaction (FSI) effects, generated 
for instance via the Bander-Silverman-Soni (BSS) mechanism J2T 



3 B —> 77 at leading power 

To leading power in AQCD/n^b and in the leading logarithmic (LL) approximation 
of QCD only one diagram contributes to the amplitude for B — > 77. It is the 
one-particle reducible (1PR) diagram with the electromagnetic penguin operator 
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b > 



s/d 




Figure 1: The leading power contribution to the B — > 77 amplitude given by the 
magnetic penguin operator Q7. The diagram with interchanged photons is not 
shown. Radiation from the 6-quark line (Fig. 0) is power suppressed. 



Q7, where the second photon is emitted from the s-quark line. An illustration is 
given in Fig. [j]. Evaluating the factorization formula (|15D for this diagram leads 
to the matrix element 



(7(*i,ei)7(*2,e2)|Q7|£) = 



(25) 



7T 



\2ki ■ e 2 k 2 ■ ei - m|ei • e 2 - 2i e(k±, k 2 , ei, e 2 



Together with 



£ 

and the definition of A± in fll2"f) this implies 
A ± = X?C 7 ^ 



X 



B 



(26) 



where CKM unitarity, + A^ = — Aj , has been used. After summing over 
the photon polarizations the branching fraction to leading power becomes 



B{B q - 77) 



T ^l92^ |A * 1 ° 7 ^Al 



(27) 



where q = d,s for the decay of a Bj or a B s meson, respectively. 

In the present approximation the strong-interaction matrix elements multi- 
plying \vp and A^ are identical and have no relative phase. The CP-violating 
quantities r^ P defined in (]20 ) therefore vanish for the strictly leading-power re- 
sult. 

One may compare the behaviour of the B — > 77 decay amplitude in the heavy- 
quark limit with that of the rare decays B — > V'-f (with V a vector meson) and 
B — > 7T7T. Disregarding CKM factors we can write 



-4(5 



77) ~ G F e 2 fBA±m 2 B 
Vj) ~Gf e mbFvm 2 B 



nil 



G F e° f n F 



B- 



m 



B 



G F e 2 A 1/2 
G F e A 3 /V- 
G F e° A 5 / 2 



5/2 



3/2 



1/2 



(28) 
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Figure 2: The subleading power 1PR diagram of the magnetic penguin operator 
Qr. The diagram with interchanged photons is not shown. 




Figure 3: The subleading power 1PI diagram. Here the four-quark operators Qi 
can contribute. The diagram with interchanged photons is not shown. 

where we have used that the form factors scale as (A = Aqcd) 

A mb f A3/2 F F^ A3/2 (W) 



We observe that the ratio of the three amplitudes in ( 28) is 1 : (A/mj) : (A/m^ 



\2 



Thus ^4(5 — > 77) is leading in the heavy-quark limit, while the other two are 
successively stronger suppressed. However, this hierarchy is compensated by the 
opposite pattern e 2 : e : 1 in the electromagnetic coupling e. 



4 Contributions to B — > 77 with power suppres- 
sion 

Subleading contributions come from the 1PR diagram, where the second photon 
is emitted off the b quark line (Fig. [J), and from the one-particle irreducible 
diagram (1PI) (see Fig. ||). 

To estimate the subleading 1PR contribution, we simply evaluate the graphs 
in Fig. fusing the 5-meson projector in (|Tfi|). The result equals Xb/ttib times the 



leading-power expression in (|25|), clearly showing the power suppression of this 
mechanism. Other corrections of the same order can arise from higher-twist terms 
in the S-meson wave function. We shall neglect all those subleading terms for 
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(b) 



1 



(c) 



lAAAA/ 




(d) 



lAAAA/ 




%0 jjtj^^^AAAA/ 




I 



(e) 



(f) 



lAAAA/ 



Figure 4: Gluon corrections to the 1PI diagram. Quark labels are as in Fig. 
and the momentum assignment is everywhere as in diagram (a). The diagrams 
with interchanged photons are not shown. 



the matrix element of Q7, keeping in mind that they could naturally contribute 
relative corrections of order ~ 10%. We stress, however, that subleading effects in 
(Qr) contribute equally to the up- and charm-quark components of the amplitude 
(see (H)). Therefore they do not give rise, in particular, to relative FSI phases 
between these sectors, which would affect direct CP violation. The same is true 
for perturbative QCD corrections to (Q7). 

We next consider more closely the diagrams in Fig. [| These 1PI contribu- 
tions, which come from the matrix elements of four-quark operators in (0), are of 
special interest for two reasons. First, they are the basic effects responsible for a 
difference between the up- and charm-quark sectors of the amplitude, including 
rescattering phases. Second, they represent the parton-level processes that are 
dual to B —>■ 77 amplitudes from Qi with hadronic intermediate states (J/^, 
D(*\ etc.), which are commonly considered as generic long-distance contribu- 
tions. Investigating the amplitude in Fig. |3] will therefore shed light on this class 
of effects in B — > 77. 

It shall now be argued that the 1PI contributions are calculable using QCD 
factorization. It is known that the one-loop contribution in Fig. |^ is free of 
infrared (IR) singularities. In addition we will show that also at 0(a s ) there 
are no collinear or soft IR divergences at leading power in l/m&. By leading 
power we here mean the lowest nonvanishing order in the power expansion, where 
the entire 1PI contribution starts only at subleading power with respect to the 
dominant mechanism in Fig. [lj This suppression will be apparent from the 
explicit expressions given in (|3(]) below. The relevant two-loop diagrams at 0(a s ) 
are shown in Fig. |j. Following the approach explained in the second reference of 
23 ] for the case of B — > Dn decays, we demonstrate the absence of divergences 
for the hard-scattering part of the amplitude at 0(a s ). As a consequence the 
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matrix elements of four-quark operators can be written in factorized form as in 
flT5|), where in the present case the kernel T^ v is independent of £: 

(l(eiHe 2 )\Qi,2\B) = fsT^e^ (30) 

The proof proceeds by identifying the potentially IR singular regions in the 
loop-momentum variables s and t and by determining the degree of IR divergence 
through power counting. We distinguish the following cases: (i) s, t soft, where 
all components scale like s, t ~ A <C m&. (ii) s, t collinear with either photon 
momentum k\ or k 2 , for instance 

s = aki + sj_ + a 2 k 2 t = (5k x + t± + (3 2 k 2 (31) 

describing the region where s and t are collinear with k\. Similar relations, with k\ 
and k 2 interchanged, apply to the case where the loop momenta become collinear 
with k 2 . The parameters in fl3~l|) scale as 



A 2 

a,/?~l, s_L,t_L~A, a 2 ,p 2 =■ (32) 

and we have 

^1,2 = ki,2 • s± = k Xt 2 ■ tj_ = (33) 

We consider explicitly the case where the quark inside the loop is massless, which 
is the most IR singular situation. The momentum I of the light external quark is 
counted as a soft quantity I ~ A. The dominant region for the diagrams in Figs. 
|3|, [| arises when both loop momenta are hard, s, t ~ mj. The diagrams then 
scale as m^. To demonstrate IR fmiteness, one has to show that the potentially 
singular regions scale at least with one power of A and thus are suppressed relative 
to the hard contribution. Scaling as hPrrib in a singular region would indicate a 
logarithmic divergence. The relevant cases for s and t are the soft-soft, soft- 
hard, hard-soft, collinear-collinear, collinear-hard, hard-collinear, soft-collinear 
and collinear-soft regions. For the collinear regions the momenta can be either 
~ k\ or ~ k 2 . In the collinear-collinear case one needs to consider primarily 
that both momenta ~ k\, or both ~ k 2 . If one momentum is ~ k\ and the 
other ~ k 2 , the contribution is always less singular in the infrared. In the second 



reference of |[24| , besides the soft-collinear regions, a supersoft-collinear scaling 
had to be discussed, where s, t ~ A 2 for supersoft momenta. In the present 
context this scaling is not more singular than the soft-collinear one and will thus 
not be considered further. 

Power counting then shows that all diagrams in Fig. |] are IR finite at leading 
power in the soft-soft, soft-hard, hard-soft, collinear-hard, hard-collinear, soft- 
collinear and collinear-soft regions. In some cases one encounters a superficial 
logarithmic divergence, which vanishes when the relations k\ 2 = and k\ • e\ — 
k 2 ■ e 2 = are used. 
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On the other hand, in the collinear-collinear region (superficial) linear diver- 
gences ~ A -1 may appear. An example is diagram (f) in Fig. |4] when both s and 
t are collinear to k\. The diagram has the form 

4 *MV- ^)7 A (l-7s)(r- fa)lv Mf+ hh a (r+ h~ /)7a(1-7s) 



D U) = d s d t 



s 2 (I - s) 2 (t - k 2 ) 2 1 2 (t + h) 2 (t + ki- s) 2 



(34) 

Expanding the numerator in powers of A using fl3l"D one finds that the superficial 
linear divergence drops out. This leads to 



D 



(/) 



cFsdH 



(l - s) 2 



(35) 



with the remainder of the diagram ~ A -7 . Together with collinear phase 
space, d A sdH ~ A 8 , the integrand then appears to scale as A°mfc. In fact, since 
fi\{j/— f) ~ iribA, (I — s) 2 ~ rrifyA for s = aki + O(A), the contribution behaves 
only as ~ A. In a similar way all other contributions from the collinear- 
collinear regions can be shown to be power-suppressed. 

In addition to the graphs displayed in Fig. [| there are also diagrams in which 
the gluon connects only to the lines inside the triangular quark loop (vertex and 
self-energy corrections). Also these diagrams are found to be IR finite at leading 
power. 

Finally, in contrast, one-gluon exchange among the external lines does gener- 
ate leading IR divergences. However, these singularities match those appearing in 
the perturbative evaluation of the I?-meson decay constant fs at the same order, 
in accordance with the factorization formula (|30D . This completes our proof that 
at 0(a s ) for the 1PI contribution there are no collinear or soft infrared diver- 
gences at leading power in 1/rrib in the hard-scattering kernel. It supports our 
claim that these diagrams are calculable using QCD factorization. 

Evaluating the 1PI diagram in Fig. |3] we find the matrix elements 




fB v ot 



71 



Qlg(zp)e(k 1 ,k 2 ,e 1 ,e2 




with 



-2 + 4,2 



Here 




and {p = u, c) 



mz 



mi 



(36) 



(37) 



(3f 



(39) 



We note that the function g(z) is related to the kernel h(u, s) defined in M in 
the context of B — > Vj by g(z) = h(l, z). Real and imaginary part of g(z) are 
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displayed in Fig. The function g(z) is regular for z — > 0, g(0) = —2. Its 
value at z = 1 is g(l) = 2(tt 2 — 9)/9, and g(z) — > as z — > oo. The function 
g(z) has an imaginary part for < z < 1/4, which is maximal close to z ~ 0.1, 
the approximate physical value for the case of an internal charm quark. The 
imaginary part is due to the rescattering process B — > pp — > 77 with an on- 
shell light-quark pair pp in the intermediate state. Since the process B — > pp is 
forbidden by helicity and angular momentum conservation for a massless quark 
p, the function g(z) is real at z = (that is for p = u). The helicity suppression 
of the phase will be absent at higher order in a s . 

Including the effect of 1PI diagrams in the decay amplitudes, the quantities 
A± from (|26|) become 

A\ = -C^ (40) 

AP _ = _ C7 ^B_h Ci + NCM) (41) 

where we have decomposed A± into 

A± = E \ { p q) A p ± (42) 

p=u,c 

We see that only A_ is modified. The imaginary part of the 1PI loop diagram 
leads to a relative CP-conserving phase between the X u and the A c contribution. 
This gives a nonvanishing r~^ P while r^ P is still zero. 

The effect of penguin operators Q3, . . . , Qq in the 1PI graphs is very small (at 
the percent level for the total amplitude) due to the small size of their Wilson 
coefficients. Also, they contribute equally to the up and charm sectors. We shall 
neglect them here. Their matrix elements can be found in J301. 
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5 Long-distance contributions to B 

D — > 77 



77 and 



We would next like to comment on the issue of long-distance effects and their 
relation with our analysis based on the heavy-quark limit. Long-distance contri- 



butions were calculated for both B — > 77 
following mechanisms have been considered: 



and D 



77 



,0- The 



B and D — > Vj, with a light vector meson V, followed by a V — > 7 con- 
version with the conversion factor supplied by the vector meson dominance 
(VMD) model. This gives a long-distance contribution to the matrix ele- 
ment of operator Q 7 . Counting powers of mj in the explicit result for the 
B s — > 07 — > 77 amplitude in |L5| shows that it is suppressed with A.Q C D/ m b 
compared to the leading power result from Fig. [I]. However, this is only one 
among the possible subleading power mechanisms for (Qj) . Another comes, 
for example, from the diagram in Fig. |^. 

We can apply a similar power counting to the D — > 77 VMD amplitude of 
3~4f. In the limit of a heavy charm quark it is likewise power suppressed with 



Aqcxj/ttt-c compared to the short- distance amplitude. The limit m c 3> J^qcd 
is certainly questionable, but has nevertheless proven successful in some 
applications to charmed hadrons [35]. Again, the VMD contributions are 



just one out of several possible power corrections to the Q7 matrix element. 

Single-particle unitarity contributions where the D° mixes with a spinless 
intermediate particle M, which then decays into a photon pair, D° — > M — > 
77. As far as the quark-level topology is concerned, such a process corre- 
sponds to the 1PI contribution in Fig. |3|. Relative to this contribution the 
amplitudes in [[34] are at most of the same order in powers of the inverse 
heavy charm quark mass. 

Two-particle unitarity contributions are other hadron-level counterparts of 
the 1PI diagrams. Triangle loops with intermediate Dy> pl|j3^| and 
34 1 mesons were calculated for B s — > 77 and D — > 77, respectively. As we 
will see in the next section, the effect of the 1PI contributions in B — > 77 



on the decay amplitude is small. The authors of [32] also get a not too 
large effect in their reestimate of the result of [51], which had suggested 
a more sizable long-distance contribution. Yet, one has to keep in mind 
that the effects from intermediate D and D* mesons are only two of many 
possible hadronic intermediate states. Since the meson model calculations 
are not based on a systematic approximation, it is for instance conceivable 
that cancellations required by quark-hadron duality would be missed. This 
could then lead to an overestimate of the effect. 



With obvious replacements we can use our 1PI results in ([36]) - ( 42|) to 
estimate the quark-level analogue of the two-particle unitarity contributions 
to D — » 77 calculated in p4| . We obtain for the 1PI amplitude 
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aD-*^ 

^ipi 



Gf ol 



V2 



71 



foQliCi + NCzyVcWVu. 




9 \™cy 
r-1 



ie(e 1 ,e2,ki,k 2 ) 



^lPI 



-2.4 ■ KT^GeV- 1 [-0.2 + 0.5i] e(ei, e 2 , fci, fc 2 ) 
1.3- lO^GeV" 1 |e(e ls e 2 ,Ari,A^)| 



(43) 



where we set fjj = 200 MeV, m s = 120 MeV, rrid = and evaluated the 
Wilson coefficients at m c . The contributions of a strange quark and a down 
quark in the loop are accompanied by CKM elements that are practically 
equal in absolute value, but opposite in sign. This leads to a strong GIM 
cancellation in the square bracket. Note that 



9{z) 



-2 + (-2 In 2 z + 2vr 2 - 4vr i In z)z + 0(z 2 



(44) 



The GIM mechanism removes the constant (—2) and gives an additional 
suppression ~ l/m? c . Due to large logarithms In z and numerical factors 
the suppression is relatively mild numerically. Also for the 1PR amplitude, 
the c —>■ wy transition exhibits a much stronger GIM cancellation than the 
corresponding process b —>■ sj. In this case this has crucial consequences for 
the hierarchy of higher-order QCD terms. Leading logarithmic QCD correc- 
tions are known to enhance the c — > w~f amplitude by more than an order 
of magnitude. Including the two-loop QCD contributions increases the am- 
plitude by another two orders of magnitude [36|. Note however that the 
two-loop estimate from |36| shows a GIM cancellation pattern very similar 
to the 1PI amplitude in (f43"D. In other words, the strong hierarchy of QCD 
corrections discussed in |36l is peculiar to the c — > transition. Hence, 



we do not expect a significant change from higher-order QCD corrections 
to ([43]). As a crude estimate, we therefore compare ( |43|) with the c — > wy 
amplitude from the dominant two-loop diagrams of which lead to an 
effective coefficient \A\ = 0.0047 (entering the amplitude with a normaliza- 
tion corresponding to IV^V^CVI). We thus obtain for the "short-distance" 
amplitude (assuming = 350 MeV): 



^SD 



4 D^77 
^SD- 



V2 7T \ D 



ie{e l , e 2 , h, fc 2 ) - k x • e 2 k 2 ■ e x + -^e x • e 2 



12.8- 10" 11 GeV" 1 |£(ei,e 2 ,A;i,fc2) 



(45) 



We see that the 1PI amplitude, in this rough order-of-magnitude estimate, 
is smaller than the CP-odd part of the short-distance amplitude. 
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CKM parameters and coupling constants 


v us 


v cb 


\v ub /v cb \ 


MS 


a 


Gp 


0.22 


0.041 


0.085 ±0.025 


225 MeV 


1/137 


1.166 x lO^GeV -2 


Parameters related to the meson 


m Bd 






^B d 


5.28 GeV 


(200 ± 30) MeV 


1.55 ps 


(350 ± 150) MeV 


Parameters related to the B s meson 










5.37 GeV 


(230 ± 30) MeV 


1.49 ps 


(350 ± 150) MeV 


Quark and W-boson masses 


m b (m b ) 


m c (m b ) 


mt(m,t) 


M w 


4.2 GeV 


(1.3 ±0.2) GeV 


166 GeV 


80.4 GeV 



Table 1: Summary of input parameters. 



The authors of [p4 |, on the other hand, find a large effect on the D —>■ 77 



branching ratio from K + K~ intermediate states. Their amplitude, how- 
ever, is formally suppressed with l/m c relative to our 1PI amplitude. In 
view of our estimate above the large effect could be overrated. Only one 
intermediate state was taken into account so that important contributions 
in the duality sum might be missing. 

The decay mechanism B s — > (pip — > 07 — > 77 was estimated to give very 
small effects on the B — > 77 decay amplitudes jl5| . It is power suppressed 
compared to the leading 1PR contribution. 



6 Phenomenology 

In this section we present numerical results for the expressions derived in this 
paper. Our choice of input parameters is summarized in Table |l|. For definiteness 
we employ the two- loop form of the running coupling a s (/i) (as quoted in 



which corresponds to a s (M z ) = 0.118 for = 225 MeV. With central values 
of all input parameters, at /i = m b , and using a nominal value of the CKM 
angle 7 = 58°, we find for the branching ratios to leading logarithmic accuracy, 
evaluating ©, © - (@|): 



B{B S -> 77) = 1.23 • 10~ 6 (46) 
B(£ d ^ 77 ) = 3.11 • 10- 8 (47) 
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To display the size of power corrections from 1PI diagrams we compare A p + with 
A p _ in (H) and (0) for B s -> 77 

A* = 4.87 A\ = 4.87 (48) 
A u _ = 5.29 = 5.07 - 0.53 i (49) 

The calculable power corrections, which manifest themselves in the difference 
between A p + and A p _, are of the expected canonical size of 0(10%). The precise 
value of these power-suppressed effects depends strongly on the renormalization 
scale 11. Typically, the relative size of the 1PI terms diminishes when the scale 
is lowered, both due to an increase in C 7 , as well as due to the simultaneous 
decrease of the combination C\ + 3C 2 . 

The main errors in (^-f4^) come from the variation of the nonperturbative 
input parameters Xb and the decay constants /# d and fs s which are all poorly 
known and enter the branching ratios quadratically. The residual scale depen- 
dence is sizeable as well because we calculated at leading logarithmic accuracy 
only. The variation of the branching ratio with the scale would be less severe 
if next-to-leading QCD corrections were known. We expect the NLL corrections 
to increase the branching ratio as was the case for both the inclusive and 
exclusive |||9] b — > 37 decays. 

The effect of next-to-leading logarithmic corrections for the branching ratio 
of b — > 57 can be reproduced by choosing a low renormalization scale (fi ~ m&/2) 
in the leading logarithmic expressions. If such a low scale were also relevant for 
B — > 77 the branching ratios would read 

77)^/3 = 1.53 - 10- 6 (50) 
B(B d -+ 77)„=W2 = 4.08 -10- 8 (51) 

However, we prefer to use the nominal /1 = and to quote the standard scale 
ambiguity with rrib/2 < fi < 2rrib as a theoretical uncertainty. 

For comparison we also show the results for the case where short-distance 
QCD effects are neglected altogether. This amounts to taking \i = M\y in the 
leading-order Wilson coefficients and gives 

B(B S 77)1*=^ = 0.60 -10- 6 (52) 
B(B d -> 77)m=m w = 1.22 -10- 8 (53) 

Clearly, the leading logarithmic QCD corrections yield a substantial enhancement 
fl5|-|T7| of the branching ratios (pEB]), (^7[), similar to the case of b — > sj. 

The numerical evaluation of (^3]) gives for central values of the input param- 
eters and at ji = m h the following predictions for the CP-violating ratios r^ p : 

rcp(B s ) = 0.35% (54) 
rcp{B d ) = -9.53% (55) 

The effect is negligible for B s — ► 77 as expected. For B d — > 77 CP violation occurs 
at the level of about 10%. It is interesting to note that in B —>■ 77 a non- vanishing 
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[l-m b /2 : 

j i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i I i i i_Z 

"" 20 40 60 80 100 120 140 160 180 

y[deg] 

Figure 6: The CP-violating ratios r~^ P for decays of neutral B& and B s mesons 
to two photons as a function of the CKM angle 7, each for three values of the 
renormalization scale \x = m&/2, and 2rrib. 

CP asymmetry appears already at 0(a®). This is possible via a modified BSS 
mechanism: Instead of the usual QCD penguin, already the electroweak loop in 
Fig. |3| itself leads to a CP-conserving phase. 

The scale dependence of the CP asymmetry is rather strong, as can be seen in 
Fig. |6|. The ratio r^p also depends sensitively on fundamental CKM parameters 
such as the CKM angle 7 (Fig. ^). The extremal value of r^ P (Bd) is obtained 
for 7 = 50 deg. The sensitivity of the branching ratios and CP asymmetries to 
variations in the relevant input parameters is summarized in Table ^. As already 
mentioned, the dominant uncertainty for our prediction of branching ratios comes 
from the variation of the hadronic parameter Xb, which enters the leading-power 
amplitude. So far, our value for \ B quoted in Table [TJ is only an educated guess 
with the uncertainty taken appropriately large. It would be highly desirable to get 
a better estimate of this important parameter. We anticipate that this problem 
can be addressed soon using QCD sum rules or lattice QCD. As As parametrizes 
the first negative moment of the B meson wave function, it is a universal quantity 
that appears in many exclusive decays. An experimental determination might also 
be possible using radiative semileptonic B decays, or even a future measurement 
of B — > 77. The variation of the remaining parameters, other than Xb, changes 
the branching ratios of B d — > 77 and B s — > 77 by about ±50% and ±35%, 



-15 
-17.5 
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B[B d -> 77) 


B(B S -> 77) 


r cp( B d) 


t Cp( B s) 




[10 8 ] 


r-i r\ fin 

[10 6 ] 


[%] 


[%] 


central 


3.11 


1.23 


-9.53 


0.35 


As 


+6.41/-1.58 


2.45/-0.61 


+4.04/-3.91 


+0.14/-0.15 


Sb 


+1.00/-0.86 


+0.34/-0.30 






IX E [rrib/2, 2m b ] 


+0.97/-0.66 


+0.31/-0.21 


+8.91/-8.07 


+0.30/-0.32 


m c 


+0.10/-0.12 


+0.03/-0.04 


+1.75/-1.51 


+0.05/-0.06 


7 = (58 ± 24)° 


+1.29/-0.96 


+0.01/-0.02 


+1.64/-0.25 


+0.06/-0.12 



Table 2: Predictions for branching ratios and CP asymmetries with the errors 
from the individual input uncertainties. 

respectively. If both branching fractions could be measured, one could consider 
their ratio where the dominant uncertainties from A# and fs would cancel to a 
large extent. The main uncertainty for the CP asymmetries at present comes from 
the variation of the renormalization scale /i. This could be reduced, in principle, 
by including higher order QCD corrections. 

Our predictions for the branching ratios are roughly two orders of magnitude 
below the current experimental bounds in (|J). The recent upper limit for B[B d — > 
77) from BaBar in ([!]) improved the previous limit from the L3 collaboration 
by a factor of 20. With more integrated luminosity accumulated the upper bound 
will be further reduced, but it will barely reach the interesting region. The high- 
luminosity option SuperBaBar suggests a total integrated luminosity of 10 ab -1 . 
The number of observed B d — > 77 events is then expected to be about two 
dozens for a branching fraction of 3 • 10 -8 [j39[] . Running such a machine at the 
T(5S) resonance would open the possibility of measuring B s —>■ 77 as well. An 
interesting option might further be the operation of a future e + e~ linear collider 
at the Z resonance. 

7 Conclusions 

In this paper we have presented a systematic discussion of the decays B s d — > 77, 
exploiting the simplifications of QCD in the heavy-quark limit. This enabled us 
to study these processes for the first time within a model-independent framework. 
In particular, we discussed the hierarchy of contributions to the decay amplitude 
in powers of Agc^/m;,. The leading contribution to B s — > 77 comes from the 
chromomagnetic b — > 57 transition where the second photon is emitted from the 
s-quark constituent. 

We have also investigated the B —>■ 77 matrix elements of four-quark oper- 
ators Qi and Q2, which proceed through 1PI triangular up- and charm-quark 
loops. These are suppressed by Aq C D/ m b relative to the leading contribution. 
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Disregarding other terms of the same order, the consideration of these 1PI dia- 
grams is nevertheless meaningful because they are the leading effects that induce 
a difference between the up- and charm-quark sectors of the amplitude. This is 
important for CP violation. They are also interesting as the partonic counter- 
parts of certain, so-called long-distance effects in B — > 77 or, similarly, D — > 77. 
We have argued that the 1PI matrix elements are calculable in QCD factorization 
and have explicitly shown this property at leading (0(a®)) and next-to-leading 
order (0(a s )) in the strong coupling. Numerically the 1PI effect is only of order 
10% in the amplitude, in accordance with a power-correction of canonical size. 

Our estimates for branching ratios and CP asymmetries are given in Table |||. 
At present there are still very large uncertainties in the absolute predictions for 
the branching fractions due to the high sensitivity on the hadronic parameter 
Xb, which is a universal, process-independent property of the B meson. The 
quantity A# is poorly studied at present, but an improved determination should 
be possible in the future. Also, a measurement of B s — > 77 could be a useful test 
of the hadronic dynamics and for A#, which plays an important role for many 
other exclusive B decays. 

The experimental study of B s ^ — > 77 will not be easy. Presumably second 
generation 5-factories or the high-intensity operation of a future e + e~ linear 
collider at the Z pole would offer the best prospects. The radiative modes B s ^ — > 
77 are interesting probes of heavy-flavour physics in their own right and especially 
in the context of factorization and the theory of exclusive B decays. In view of 
this, new ways towards their experimental observation should continue to be 
pursued. 
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